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ABSTRACT 


Frequency domain parameter identification techniques provide a straightforward 
approach to transfer function estimation. However, for high-order systems, numerical 
difficulties may be encountered during the estimation process. Inaccuracies may result 
because of the large variation of the transfer function polynomial coefficients for high- 
order systems. The lack of numerical precision to represent this variation may cause 
the estimation process to break down. 

This paper presents a technique for estimating transfer functions in partial fraction 
expansion form from frequency response data for a high-order system. The problem 
formulation avoids many of the numerical difficulties associated with high-order 
polynomials and has the advantage of having the option to fix the damping and 
frequency of a mode. If known, during the estimation process. The resulting transfer 
funclion(s) may be converted to Jordan-Form time domain equations directly. 

During the implementation of this technique, a frequency and amplitude normalizing 
window was developed that maximized the efficiency of the optimization algorithm. 
The combination of estimating the transfer function In factored form, the ability to fix 
previously determined parameters and the effectiveness of the normalizing window led 
to a progressive approach to synthesizing transfer functions from frequency response 
data for high-order systems. 
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Frequency domain parameter identification techniques provide a straightforward 
approach to transfer function estimation. However, for high order systems, 
numerical difficulties may be encountered during the estimation process. 
Inaccuracies may result because of the large variation of the transfer function 
polynomial coefficients for high order systems. The lack of numerical precision 
to represent this variation may cause the estimation process to break down. 

This paper presents a technique for estimating transfer functions in partial 
fraction expansion form from frequency response data for a high order system. 
The problem formulation avoids many of the numerical difficulties associated 
with high order polynomials and has the advantage of having the option 
to fix the damping and frequency of a mode, if known, during the estimation 
process. The resulting transfer function(s) may be converted to Jordan-Form 
time domain equations directly. 

During the implementation of this technique , a frequency and amplitude normal izing 
window was developed that maximized the efficiency of the optimization algorithm. 
The combination of estimating the transfer function in factored form, the 
ability to fix previously determined parameters and the effectiveness of 
the normalizing window led to a progressive approach to synthesizing transfer 
functions from frequency response data for high order systems. 
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NUMERICALLY EFFICIENT ALGORITHM FOR 
MODEL DEVELOPMENT OF HIGH ORDER SYSTEMS 
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cycle 

Statistical methods applicable 


PREVIOUS WORK 
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Transformation errors 

Sensitivity of polynomial roots to variations in polynomial coefficients 


SUMMARY OF CONTRIBUTIONS 
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Stepwise approach for synthesizing transfer functions where order of sys- 
tem is high and unknown 


FACTORED FORM ESTIMATION 
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where: N = order of system 

Q = # of second order terms 



Express measured and estimated frequency responses in terms of real and imaginary 
components 
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Solve for unknown parameters 

Set partial derivatives equal to zero 

Solve using nonlinear optimization technique 


Fifth Order Single Precision Example 

Simulate parameter identification of high order system 
Modes distributed over wide frequency range 
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DENOMINATOR COEFFICIENTS 
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= Single Precision Variable Representation 
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Iterations converge to minimization of | e 





FACTORED FORM APPROACH 
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Benefits - Need for strong initial guess eliminated 

effectiveness of optimization algorithm maximized 



STEPWISE FACTORED FORM TECHNIQUE 



Scale Amplitude Data 
around 1 (abs.) 


Identify 

Parameters 


Improve x 
Parameters, 
\ ? / 


Parameter = 

Estimated 

Vector 

Values 


.1 = 1+1 



Run Nth Order Case > 
over entire frequency 
range using estimated 
values as initial 
Parameter Vector J 








5th Order Example: Fi 
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SCALING: 

FREQ = 1.0 
AMP = 1 X 10" 3 
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CONCLUSIONS 
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Development of frequency and amplitude normalizing window that eliminates the 
initial guess problem and maximizes the effectiveness of the optimization algo- 
rithm 


